Deeply Virtual Compton Scattering off the Neutron 
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The present experiment exploits the interference between the Deeply Virtual Compton Scattering 
(DVCS) and the Bethe-Heitler processes to extract the imaginary part of DVCS amplitudes on 
the neutron and on the deuteron from the helicity-dependent D(e, e'y)X cross section measured at 
Q 2 =1.9 GeV 2 and xb=0.36. We extract a linear combination of generalized parton distributions 
(GPDs) particularly sensitive to E q , the least constrained GPD. A model dependent constraint on 
the contribution of the up and down quarks to the nucleon spin is deduced. 

PACS numbers: 13.60.Fz, 13.85.Hd, 14.20.Dh, 14.65.-q 



Understanding the structure of the nucleon in terms of 
quarks and gluons is a central project of modern hadronic 
physics. In the non-perturbative regime relevant to nu- 
clear scales, Quantum Chromodynamics (QCD), the the- 
ory describing the elementary dynamics of the nucleon, 
is not yet solvable and remains rather mysterious. The 
electromagnetic probe provides an outstanding tool to 
study the nucleon structure. In this letter, we present 
the first study of the (e, e'7) reaction on neutrons off a 
deuterium target. 

Elastic electron scattering revealed the non-pointlikc 
nature of the nucleon 111 . Deviations from the Mott cross 



section define the electromagnetic form factors which de- 
scribe the spatial distribution of charge and current in- 
side the nucleon, as functions of the invariant momentum 
transfer squared (Q 2 ), i.e. the resolution of the probe. 
Deep Inelastic Scattering (DIS) revealed the partons in- 
side the nucleon The DIS cross section can be ex- 
pressed in terms of the probability to find a quark with 
fraction xb of the nucleon longitudinal momentum. This 
motivated extensive measurements of the momentum dis- 
tribution of quarks and gluons in nucleons, i.e. the parton 
distributions. Polarized DIS with longitudinally polar- 
ized beams and targets measures the probability to find 



a parton of given momentum with spin aligned or anti- 
aligned with the proton spin. These experiments yielded 
the unexpected result that the quarks carried about only 
30% of the total spin of the nucleon [§| and questioned 
the role of gluons in this puzzle. Experimental results 
to date suggest that the gluon polarization does not con- 
tribute significantly to the nucleon spin 0, @] • A natural 
candidate to solve this problem is the orbital angular mo- 
mentum of quarks and gluons. 

The Generalized Parton Distribution (GPD) frame- 
work provides a new formalism to unravel the nucleon 
structure that unifies form factors, parton distributions, 
and the angular momentum of partons 0, 0, @| ■ Quark 
GPDs are four universal functions H q , E q , H q , and E q , 
defined by the nucleon hclicity-conserving and helicity- 
flip matrix elements of the vector and axial-vector cur- 
rents for quark flavor q. The GPDs correspond to the 
amplitude for removing a parton of momentum fraction 
x + £ and restoring it with momentum fraction x — £ 
(Fig. [I]). In this process, the nucleon receives an invariant 
momentum transfer t—A 2 . In impact parameter space, 
the transverse momentum transfer is Fourier conju- 
gate to the transverse position of the parton leading to a 
femto-tomography 0, 1(J EH E3 of the nucleon: GPDs 
represent distributions, in the transverse plane, of par- 
tons carrying longitudinal momentum x. The correlation 
of the position and momentum of quarks can exactly be 
combined into an orbital momentum [131 ] . This is explicit 
in Ji's sum rule [8j 



Jq — 2^^9 "F 



(1) 



1 f +1 

= - J ^ dxx[H q (x,£,t = 0) + E q (x,Z,t = 0)] 

which is independent of the £-value. Polarized DIS mea- 
sures the spin part AE 9 , and unpolarized DIS determines 
the momentum sum rule and forward limit of H a 



dxxq(x) = / dx xH q (x,£ = 0,t = 0) (2) 



M q = 



meaning that the contribution of H q to Eq. [T] is known. 
Constraints on E q will allow access to the quark orbital 
momentum in the nucleon. 

Deeply virtual Compton scattering (DVCS) is the sim- 
plest reaction to access GPDs (Fig. [T]). In the Bjorken 
limit, similar to DIS, where —t <C Q 2 and Q 2 is much 
larger than the quark confinement scale, the factoriza- 
tion theorem separates the reaction amplitude into the 
convolution of a known perturbative j*q — > 79 kernel 
with an unknown soft matrix element describing the nu- 
cleon structure (GPDs) The Bethe-Heitler (BH) 
process, where the real photon is emitted by either the 
incoming or scattered electrons, serves as a reference am- 
plitude that interferes with the Compton amplitude. The 
difference between polarized cross sections for opposite 
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FIG. 1: Lowest order (QCD) amplitude for the virtual Comp- 
ton process. The momentum four-vectors of the incident and 
scattered photon are q and q' , respectively. The momentum 
four- vectors of the initial and final proton are p and p , with 
A=(p'-p)=(q-q'). The DIS scaling variable is x B =Q 2 /(2p-q) 
and the DVCS scaling variable is £=:eb/(2 — xb). In light cone 
coordinates defined by P=(p + p')/2, the initial and final mo- 
mentum of the photons are — 2£ and 0, respectively. 



beam longitudinal polarization isolates, at leading order 
in 1/Q, the imaginary part of the interference between 
the BH and DVCS amplitudes This difference is a 



direct measurement of a linear combination of GPDs [17 1 
dominated by the contribution of E q in the neutron case. 

The first evidence for DVCS was reported in beam- 
helicity asymmetries at HERMES [Hj], and CLAS 18 1, 
and in unpolarized cross sections at HERA 2(| 21, 22 j, 



and more recently in measurements of the tar get spin 
asymmetry [23| and beam charge asymmetry 24 [. A 
dedicated H(e, e'j)p experiment in Jefferson Laboratory 
(JLab) Hall A shows evidence for factorization at Q 2 as 
low as 2.0 GeV 2 j2j|. The E03-106 experiment [26| re- 
ported here is an exploratory experiment investigating 
for the first time the DVCS reaction off the neutron. 

The experimental data were acquired in JLab Hall 
A, consecutively to the H(e, e'j)p experiment. A 
5.75 GeV/c longitudinally polarized electron beam im- 
pinged on a 15 cm liquid D2 cell serving as quasi-free neu- 
tron target. Scattered electrons of 2.94 GeV/c were de- 
tected at 19.3° in the left High Resolution Spectrometer 
(HRS-L) [13 selecting kinematics at Q 2 =1.9 GeV 2 and 
xb=0.36. DVCS photons were detected in a PbF2 elec- 
tromagnetic calorimeter organized in an 11x12 array of 
3x3x18.6 cm 3 crystals centered around the direction of 
the virtual photon at —18.3°. The calorimeter front face 
was 110 cm from the target center covering a t acceptance 
—0.5 GeV 2 < t. Typical beam intensities of 4 ^A yielded 
a 4xl0 37 cm _2 -s _1 /nucleon luminosity with 76% polar- 
ized electrons. Three independent reactions are used to 
calibrate and monitor the calorimeter: H(e, en,i„ Phrs); 

D ( e > e Calo. 7r HRs)PP> andH > D ( e . e HRS 7r Calo.) X O- Recoil 

particle detection in the HRS provides tagged electrons 
in the calorimeter allowing for the independent calibra- 
tion of each block. The mass and width of the ir° peak 
reconstructed from the invariant mass of 77 events in the 
calorimeter provide independent tests of the previous cal- 
ibrations. 7Tjj Rg and 7TQ a | data have been taken simulta- 
neously with DVCS data, ensuring a continuous monitor- 
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FIG. 2: (top) Helicity signal (Eq. 0} for D(e,e'7)X and 
H(e, e /r y)X events; H2 data are folded with a momentum dis- 
tribution of the proton in deuterium, and scaled to the D2 
data luminosity; the simulation curve is for the Fermi broad- 
ened H(e, e'~/)p reaction, (bottom) Residual helicity signal 
after H2 subtraction; the arrows indicate the M\ average po- 
sition of n-DVCS and d-DVCS events for < t >=-0.3 GeV 2 ; 
the simulation curves, integrated over the complete exper- 
imental acceptance and obtained for the arbitrary values 
3m [Cn] exp =-Qm [C^^^-l (Eq. HJ, illustrate the sensi- 
tivity of the data to the neutron and deuteron signals. 



ing of the calibration and the resolution of the calorime- 
ter. A 1% uncertainty on the calorimeter calibration was 
estimated from the differences between ir~ and ir° cali- 
brations. The final state of the D(e, e'j)X reaction was 
selected via the squared missing mass M^ = (q + p — q') 2 
reconstructed from the virtual and real photons. 

The three- momentum transfer |A| to the target varies 
within 0.4-0.8 GeV/c in our acceptance. In this range, 
the impulse approximation (IA) is expected to accurately 
describe the inclusive yield. Within the IA, the cross sec- 
tion for electroproduction of photons on a deuterium tar- 
get may be decomposed into elastic (d-DVCS) and quasi- 
elastic (p-DVCS and n-DVCS) contributions following 

D(e, e'^)X = d(e, e'j)d+n(e, e'j)n+p(e, e r f)p+. . . (3) 

where meson production channels are also contributing as 
background. Cross sections are obtained from D(e, e'j)X 
events after subtraction of the proton quasi-elastic contri- 
bution deduced from measurements on a liquid H2 target: 
the Fermi motion of bound protons is statistically added 
to the squared missing mass -M^-jo of free proton data 
following M\—M\ |o — 2pi ■ (q — q 1 ) where fj is the initial 
proton momentum in the deuteron from [291 ] ; this leads 
to a 3% relative increase of the M\ spectrum resolution. 
The helicity signal (Sh) is defined according to 



S h = 



(AT H 



)d 5 $- f 



(N + -N-)d 5 <I> (4) 



where d 5 $ = dQ 2 dx-Bdtd4> e d4> 11 is the detection hyper- 
volume; the integration boundaries in Eq. [4] define the 
limits in the azimuthal angle </> 77 [30j; are the num- 
ber of counts for ± beam helicity, corrected for random 
coincidences, and integrated over a particular bin in Mj^ . 
The helicity signal for D2 and H2 targets from (e, e'j) co- 
incident detection is displayed in Fig. [2] (top) as a func- 
tion of the squared missing mass. For our purposes, M\ 
is calculated with a target corresponding to a nucleon 
at rest, leading to the kinematic AMjf ~ t/2 separa- 
tion between deuteron elastic and nucleon quasi-elastic 
contributions. Pion production channels (eA — > eA^-n, 
eA — ► eATr°Tr . . .) are strongly suppressed by the kine- 
matical constraint M\ < (M + m w ) 2 =Mj (: \ cnt . Their 
contribution to the helicity signal of p-DVCS, induced 
via resolution effects below M^-| cut , was found to be neg- 
ligible on the proton as illustrated by the comparison 
between H2 data and scaled simulations (Fig. [2] top). 
Figure [2] (bottom) shows the subtraction (D— H data) 
of the two spectra of Fig. [2] (top). The residual helic- 
ity signal for M\ < Mx\ C ut is compatible with zero. It 
corresponds to the sum of the coherent d-DVCS and in- 
coherent n-DVCS processes (Eq. [3]) . Asymmetric decays 
of 7T° (in eA — > eAir ), where only one photon is detected 
in the calorimeter, mimic DVCS events. The contamina- 
tion due to this background was treated as a systematic 
error estimated from the number of detected ir° events, 
corresponding to primarily symmetric decays [2F3 | . 

The H2 results [25( show that the handbag mechanism 
(Fig. QJ dominates the p-DVCS helicity-dependent cross 
section difference at our kinematics. As a consequence, 
only twist-2 contributions are considered in this analy- 
sis. The exp superscript in Eq. [5] reflects this restriction. 
In the impulse approximation, we write the experimen- 
tal helicity-dependent cross-section difference as the sum 
of the (incoherent) neutron and the (coherent) deuteron 
contributions, within the formalism of Refs. [l?], l3ll | 
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d are kinematical factors with a 



77 dependence that 
arises from the electron propagators of the BH amplitude; 
in the t range of interest, the averaged /T d ratio vary 
from 0.4 to 0.9 with increasing \t\. 3m [C^] depends on 
the interference of the BH amplitude with the set T = 
{H, £, H} of twist-2 Compton form factors (CFFs): 



[c I n } exp ^[c I n ] = F 1 H + m + F2)n 



t 



AM 2 



F 2 £ (6) 



where -F^i 7 ^) is the Dirac(Pauli) form factor entering into 
the BH amplitude. Similarly, 3m [C d ] depends on the 
different set of spin-1 CFFs of the deuteron 31]. The 
imaginary part of twist-2 CFFs is determined by the x = 
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FIG. 3: The t-dependence of the extracted sin(<^ 77 ) moments 
for coherent d-DVCS (top panel) and incoherent n-DVCS 
(bottom panel). Error bars show statistical uncertainties; sys- 
tematical uncertainties are indicated by the shaded bands. 



±£ points of the GPDs, with for example: 

3m [£] = rr ]T e 2 (E q (t £, t) - E q (-£, £ t)) 



(J) 



where e q is the quark charge in units of the elementary 
charge. While Eq. [6] for a proton is dominated by H and 
TL 1 it becomes essentially sensitive to £ in the neutron 
case following the small value of F\ and the cancellation 
between u and d polarized parton distributions in 7i [32| . 

3m [C^] exp and 3m [C^] exp are simultaneously ex- 
tracted in each i-bin from a global analysis involving 
7x12x30 bins in t <g> 77 <g> M| 6 [-0.5; -0.1] GeV 2 <g> 
[0;2tt] <g) [0.;1.15] GeV 2 . A Monte Carlo simulation 
with the kinematic weights of Eq. [5] as a function of 
(t, <^ 7 -y, Af|) is fitted to the experimental distribution 
[iV+(t, 7T ,M|-) -N~(t, </>iy,M%)] obtained after the 
D— H subtraction. The two coefficients 3m [C^(tj)] exp 
and 3m [C^(t$)l exp are the free parameters of the fit in 
each bin t i . The binning in 77 allows the determination 
of the sin(</> 77 ) moments whereas the binning in M\ al- 
lows the separation of the d-DVCS and n-DVCS signals. 
The simulation includes both external and real internal 
radiative effects. It takes also into account detector res- 
olution and acceptance. Finally, virtual and soft real 
radiative corrections are applied with a global correction 
factor of 0.91±0.02 to the experimental yields [33| . 

Figure [3J displays the experimental values (Tab. [J) of 



exp 



At low |£|, the small kinematic sepa- 
ration between d-DVCS and n-DVCS leads to a strong 
anti-correlation between deuteron and neutron moments 
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-0.69 ± 0.38 ± 0.24 


0.98 ± 0.57 ± 0.33 
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-0.225 


0.67 ± 0.48 ± 0.39 


-1.22 ± 0.69 ± 0.40 


-0.91 


-0.166 


-1.54 ± 0.80 ± 0.52 


2.32 ± 1.04 ± 0.61 


-0.95 



TABLE I: Experimental values of the sin(</> 77 ) moments as a 
function of t (in GeV 2 ). The first error is statistical and the 
second is the total systematic one resulting from the quadratic 
sum of each contribution; a n d is the correlation coefficient 
between the two extracted moments. 



(Tab.Q}. The larger statistical errors on the extraction at 
low \t\, in spite of higher absolute statistics, reflect this 
feature. The systematical errors come essentially from 
the t-dependent uncertainties on the relative calibration 
between D2 and H2 data, and estimates of the bound 
on 7T° contamination; other contributions originate from 
DVCS detectors acceptance and luminosity (3%), beam 
polarization (2%) and radiative corrections (2%). As ex- 
pected from Fig. the moments are globally compatible 
with zero. Experimental results are compared to model 

37 1 



34, uiol and neutron [3i 



calculations for deuteron 
GPDs. The deuteron calculations exhibit a rapid de- 
crease of the deuteron form factors with \t\. The n-DVCS 
results are compared to two different models: one where 
the GPDs parametrization is constrained by lattice cal- 
culation of GPDs moments 36], and another where E q 
is parametrized by the unknown contribution of valence 
quarks to the nucleon angular momentum (32| . Both ap- 
proaches reproduce the rather flat i-dependence of the 
data. Three examples of calculations corresponding to 
different values of the u (J u ) and d (J<j) quark contri- 
butions are shown. This comparison indicates that the 
present data provide constraints of the GPD models, par- 
ticularly on E q . 

A correlated constraint on J„ and Jd can be extracted 
from a fit of the VGG model 31 
(Fig. [3]) , relying on the y 2 quantity 



37j to the neutron data 
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(3m 



(Ssys) 
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The condition X 2 < Xmm + 1 (Xmm/ D oF=6.6/5) defines 
the band J d + (J J 5.0) = 0.18 ± 0.14 of Fig. H The 
model dependence of this analysis should be stressed: n- 
DVCS data involve GPDs at one point £=±£ and t ^ 
while the Ji sum rule (Eq. [TJ) is an integral over x ex- 
trapolated to t=0. A similar constraint obtained from 
HERMES preliminary p-DVCS data on a transversely 
polarized target 38|, |39( is also shown together with lat- 
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FIG. 4: Experimental constraint on J u and Jd quark angular 
momenta from the present n-DVCS results. A similar con- 
straint from the p-DVCS target spin asymmetry measured by 
HERMES [31J391, and different lattice QCD based calcula- 
tions 36, 40, 41, 42] arc also shown. 



tice QCD based predictions SI, EEL 51. it remains a 
future theoretical study to evaluate, in a model indepen- 
dent way, the constraints on J u and Jd from a finite set 
of measurements. As expected from isospin symmetry, 
n-DVCS data have enhanced sensitivity to the d quark 
of the proton relative to p-DVCS data. This complemen- 
tarity is a key feature for future experimental programs 
investigating quark angular momenta. 

In summary, this experiment provides a determination 
of the i-dependence of a linear combination of GPDs 
from the n-DVCS helicity-dependent cross-section differ- 
ence in the range [—0.5; —0.1] GeV 2 . These data, mostly 
sensitive to E q , were found to be compatible with zero. 
The coherent d-DVCS contribution was also extracted, 
the high \t\ behaviour being compatible with expecta- 
tions. We provide the first experimental constraint on the 
parametrization of the GPD E q that can be expressed, 
within a particular model, in terms of a constraint on 
the quark angular momenta. DVCS experiments on the 
neutron appear as a mandatory step towards a better 
knowledge of the partonic structure of the nucleon. 
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